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ABSTRACT 


A  design  method  for  solving  the  problem  of  robustness  to 
cross-coupling  perturbations  in  multi variable  control 
systems   for      the    X22A     V/SIOL    aircraft      is    presented.  The 

method  uses  numerical  optimization  procedures  to  manipulate 
the  system  feedback  gains  as  direct  design  variables.  The 
manipulation  is  accomplished  in  a  manner  that  produces 
desired  performance  by  pole  placement  and  robustness  by 
modification  of  the  minimum  singular  values  of  the  system 
return    difference    matrix. 

Channels  affected  by  cross-coupling  perturbation  may  be 
recognized  by  the  character  of  their  transfer  function  Bode 
plots.  The  mechanism  used  by  the  pole  placement  and  robust- 
ness routine  in  obtaining  a  robust  design  is  evident  from 
the  gain  changes  associated  with  the  transfer  function 
diagram  and  the  zero  shifts  shewn  on  pole-zero  plots.  The 
pole  placement  and  robustness  routine  uses  gain  equalization 
and  zero  assignment  to  modify  the  characteristics  of  the 
system  in  the  areas  of  low  singular  values,  producing  a 
robust   design. 
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I.  INTRODUCTION 

In  practice,  the  control  system  designer  carely  has  the 
freedom  to  feed  back  the  entire  state  of  the  system  to 
achieve  the  desired  system  performance.  Thus,  an  important 
problem  from  a  practical  standpoint  is  the  determination  of 
constant,  output  feed  back  gains  for  the  control  of  systems 
with  unmeasured  states. 

The  pole  placement  only  and  pole  placement  and  robust- 
ness design  procedure  uses  placement  to  establish  a  designer 
selected  performance  level  and  then  a  minimum  singular  value 
level  to  establish  robustness.  The  pole  placement  only  was 
relatively  easy  to  implement  through  a  numerical  optimiza- 
tion routine.  By  using  this  numerical  procedure  it  is  also 
simple  to  incorporate  robustness  into  the  procedure  along 
with  the  performance  requirement.  The  technique,  which 
utilizes  a  modern  optimization  routine,  can  significantly 
assist  the  designer  in  obtaining  robustness  in  the  face  of 
cross  coupling  perturbations.  It  has  been  shown  that  the 
cross  coupling  problem  can  be  detected  by  using  classical 
open  loop  Bode  diagrams  as  well  as  modern  control  analysis 
based  on  singular  values.  As  currently  employed,  the  pole 
placement  and  robustness  program  is  used  to  obtain  pole 
placement  and  robustness  for  a  given  set  of  starting  gains 
and  a  selected  optimization  routine  from  the  ADS  program. 

The  pole  placement  and  robustness  design  procedure  is  a 
straight  forward  numerical  optimization  procedure  for  the 
practical  application  modern  MIMO  system  analysis.  The  new 
aspects  of  the  procedure  are  the  implementation  of  both  pole 
placement  and  robustness  criteria  within  the  same  design 
program. 
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The  pole  placement  and  robustness  design  routine  devel- 
oped for  this  thesis  has  been  used  on  X-22A  A/C  problem.  In 
these  studies  the  pole  placement  and  robustness  design  code 
has  proven  capable  of  meeting  the  desired  goals  of  pole 
placement  and  robustness  and  also  brought  to  light  some 
interesting  aspects  of  the  cross-coupling  perturbation 
problem. 

The  remainder  of  the  thesis  will  present  background 
material  on  Robustness (SISO  system)  in  Chapter  Two  and 
Robustness  of  MIMO  systems  in  chapter  Three.  Pole  placement 
only  design  procedure  in  chapter  Four,  Pole  placement  and 
Robustness  design  in  chapter  Five,  along  with  a  discussion 
of  Applications  for  X-22A  a/c  problem  by  input  transfer 
function  F*G  analysis  and  output  transfer  function  G*F  anal- 
ysis. Conclusions  will  be  presented  in  the  final  chapter. 
The  computer  programs  used  in  the  present  analysis  were 
developed  in  [Ref.  9  ]- 
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II.  ROBUSTNESS:  SIHGLE^INPOT  S INGLE-OOTPDT  SYSTEKS  (SISO) 

A  review  of  the  concept  of  robustness  and  stability  in 
the  framework  of  a  conventional  SISO  system  will  be  done 
before  pursuing  the  concepts  in  a  more  complicated  fashion 
in  the  following  chapters.  A  simple  interpretation  of 
robustness  is  the  ability  of  the  system  to  tolerate  design 
perturbations.  These  perturbations  could  be  in  the  form  of 
actuator  failures,  plant  parameter  uncertainty,  unmodeled 
dynamics  or  nonlinear  terms,  or  any  one  of  many  otuer 
perturbations  to  the  nominal  design  of  the  system. 

The  primary  reason  for  feedback,  systems  is  the  control 
of  uncertainty  within  the  system.  By  appropriate  use  of 
feedback,  properties  that  would  lead  to  an  unstable  system 
may  be  controlled.  When  stability  and  robustness  aspects 
are  considered  for  a  SISO  system,  frequency  domain  design 
concepts,  using  either  Nyguist  or  Bode  plots,  are  normally 
used.  Robustness  in  SISO  systems  is  formulated  naturally  by 
the  concept  of  gain  and  phase  margins,  both  of  which  are 
readily  available  on  the  Nyguist  or  Bode  diagram. 
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Figure    2.1  .      Classical    Bode   Plot. 

Gain  and  phase  margin  can  be  defined  in  terms  of  the 
open-loop  frequency  domain  plots  in  either  the  Bode  or 
Nyquist  format.  Figure  2.1  depicts  a  classical  Bode  plot 
showing  gain  and  phase  margin  detarmina tion  from  the  plot. 
The  Nyguist  plot  may  also  be  used  to  obtain  this  informa- 
tion. Nyquist  criterion  states  that  if  the  open-loop 
transfer  function  G(s)H(s)  has  n  poles  in  tha  right  half 
plane  and  the  limit  of  G  (s)  H  (s)  =constant  as  s— »oothen  for  a 
stable  system  the  locus  of  G(s)H(s)  will  encircle  the  -1+jO 
point  n  times  in  the  counterclockwise  direction  as  s  varies 
along  the  Nyquist  contour.  If  there  are  no  poles  in  the 
right  half  s  plane  then  the  locus  will  not  encircle  the  -1  + 
jO  point.  The  diagram  in  Figure  2.2  illustrates  a  nominally 
stable  system.  The  gain  and  phase  margin  may  be  determined 
directly   from   the   diagram.- 

Any  change  in  the  loop  transfer  function,  provided  the 
order  of  G(s)H(s)  does  not  change,  that  changes  the  number 
of    times   the   locus   of   G(s)H(s)     encircles    the    (-1/0)     point   in 
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Figure   2.2        Nyquist   Plot   of   Stable   System. 

the  Nyquist  plot  causes  the  system  to  become  unstable.  This 
leads  to  the  conclusion  that  the  minimum  distance  of  the 
locus  of  G(s)H(s)  to  the  (-1,0)  point  is  a  measure  of  the 
system  stability.  This  distance  concept  carries  over 
directly  to  the  Multi-Input  Multi-Output  (MIMD)  system  as 
will  be  shown  in  the  next  chapter.  Examples  of  a  multipli- 
cative perturbation  and  an  additive  perturbation  illustrate 
this  idea.  Figure  2.3  is  an  additively  perturbed  system. 
Figure  2.4  shows  the  Nyquist  plot  for  this  system.  Assuming 
that  the  plant  is  itself  stable  and  the  perturbations  are 
also    stable   the    diagram   may    then  be    used    to    determine   how 
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Figure  2.3    Additively  Perturbed  System. 

near  the   system  is  to  instability  for  the   given  perturba- 
tion. 

Since  the  system  is  stable  the  (-1/0)  point  is  encircled 
the  correct  number  of  times  by  the  nominal  plant.  If  the 
locus  of  g(jCW)  in  the  diagram  is  warped  until  it  passes 
beyond  the  (-1,0)  point  then  clearly  the  number  of  encircle- 
ments of  this  point  will  change  and  the  system  will  become 
unstable,  assuming  the  order  of  the  plant  is  not  changed  by 
the  perturbation.  To  keep  the  locus  of  points  from  moving 
beyond  the  (-1,0)  point  equation  2.1   must  hold. 
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Figure    2.4        Additi-ve    Nycjuist   Plot, 


I  Ag(jco)  l<l    1    +   g  (jcu)  | 


(2.1) 


This  condition  is  illustrated  in  Figure  2.5.  The  right-hand 
side  of  equation  2.  1  is  just  the  magnitude  of  the  return 
difference  transfer  function  of  the  nominal  system.  The 
multiplicative  case  is  depicted  in  Figure  2.6  with  its  asso- 
ciated Nyguist  plot  in  Figure  2.7.  The  requirement  for 
stability  is  similar  to  the  additive  case  and  may  be  stated 
in  equation  2.2 


1  Ag(jco)  I  <l  1  ♦  (g(jw)  -1)  I 


(2.2) 
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Figure  2-5    Hyquist  for  Inequality  Additive  Condition. 

The  above  arguments  will  be  applied  again  in  Chapter  3  to 
develop  multivariable  stability  and  robustness  properties. 
With  this  basic  review  of  the  concepts  of  stability  and 
robustness  in  the  classical  SISO  system  complete,  the  next 
chapter  will  extend  some  of  these  basic  concepts  to  the  dIMO 
system. 
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Figure  2.6   Multiplicative  System. 
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Figure  2.7    Nyguist  Plot  for  Multiplicative  Systea, 
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III.  HULTIVARIABLE  SYSTEMS  ROBUSTNESS 

A  generalization  of  the  SISC  Hyguist  theory  discussed  in 
the  previous  chapter  has  been  made  for  the  HIHG  prollem. 
This  generalization  leads  directly  to  the  application  of 
singular  value  concept.  The  generalization  is  expressed  in 
the  form  of  the  muitivariable  Nyquist  theorem  which  requires 
that  a  closed  loop  stable  system  have  the  same  number  of 
counterclockwise  encirclements  of  the  origin  by  the  locus  of 
the  det  (I+G  ( ju>)  )  as  the  number  of  open  loop  poles  that  are 
in  the  right  half  plane.  This  theorem  is  formally  stated  as; 
let  N[f  (s)  ]  denote  the  number  of  clockwise  encirclements  of 
(-1,0)  by  the  locus  of  f  (s)  as  s  traverses  the  contour  D 
of  Figure  3.1  in  a  clockwise  sense.  The  closed-loop  system 
will  be  stable  if  and  only  if  for  all  R  sufficiently  large 

N[f  (s)  ]  =  -P 
where      N  =  number  of  encirclements 

f  (s)  =  -1  +  det[I+G  (s)  ]  =   f£_  (s)/  (pcu    (s)  -1  and 

P   =  the   number   of  closed   right-half   plane  zeros  of 

9.L  <s>- 

The  application  of  the  Nyquist  theorem  comes  through  the 
fact  that  a  muitivariable  system  will  not  be  robust  to 
modelling  errors  if  the  return  difference  matrix,  I+G,  is 
nearly  singular  for  some  frequency.  If  I  +  G  is  nearly 
singular  a  small  change  in  G  may  make  I  +  G  exactly 
singular.  This  causes  the  det (I  +  G)  to  become  zero  and  the 
Nyguist  encirclement  count  to  change  indicating  an  unstable 
system.  It  is  possible  for  very  small  changes  in  I+G  to 
produce  large  changes  in  the  determinantof  1+3.  The  matrix 
I+G 
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Figure   3-1         Nyguist   D    Contour. 


10/s+a  9.99/s+a 

10/s+a  10/s+a 

has  determinant  0.1/ (s+a)2.  If  the  element  p  is  changed 
by  only  one  percent  to  9.9/s+a  the  determinant  becomes 
1.1/ (s+a)2  which  is  a  significant  change  in  the  determinant- 
value.  Therefore,  it  is  evident  that  det  (I  +  G)  is  not  an 
accurate  measure  of  how  near  the  return  difference  is  to 
singularity.  Researchers,         in      the        field      of      controls 

[Ref.  1],  [Bef.  2],  [Hef.  3]#  [fief.  4]  have  used  singular 
value  analysis  to  determine  how  near  the  return  difference 
matrix   is  to   singularity. 

Since   the   number    of   encirclements  of   the   Nyguist    diagram 
changes   as  f (s)     passes   through    the   -1    point    or    when    det(I+G) 
is    zero   it  is   important   to      find  how    near    the  return    differ- 
ence   matrix      I+G    is      to   being      singular.         This      nearness   to 
singularity   can    be    interpreted    as    the   distance    of    matrix    I+G 
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to  the  critical  point,  -1.  A  guantity  which  can  be  used  to 
express  the  nearness  to  singularity  of  the  matrix  is  the 
minimum  matrix  singular  value  denoted  by^s  .  Given  a  matrix 
A   the    singular    value    may    be    expressed   by    eguation    3. 1 


0\    (A)    =    min    (    A;.  (U)) 


(3.1) 


where  J\-   (A  A)   is   the  eigenvalue  of  the 
transpose  of   A,  times  A.    A   basic  MIHO 


complex  conjugate 
linear  system   is 


Figure  3.2    Basic  Multi-input   Multi-oatput  System. 

depicted  in  Figure  3.2.  An  additive  perturbation  to  the 
plant  is  shown  in  Figure  3.3.  If  the  plant  is  stable  before 
the  perturbation  is  added  to  the  system  the  Nyguist  theorem 
will  be  satisfied  and  the  locus  of  GH  will  not  encircle  the 
-1,0  critical  point.  When  the  perturbation  is  added  to  the 
system   as   long  as   the   Nyguist   locus   is  not   forced  to 
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Figure   3.3        Additive   Perturbation. 

encircle  the  -1  point  the  system,  will  remain  stable.  A 
sufficient  condition,  recalling  the  SISO  discussion  in 
chapter  2,  for  the  perturbed  Nyguist  plot  not  to  change 
encirclements  is  that  the  norm  of  the  perturbation  Ag 
remain  less  than  the  norm  of  the  return  difference  matrix  as 
expressed   in   equation   3.2. 


II    A£(jOJ)||<    1/||  (I+G)  -i)  11  (3.2) 

The    condition     uj  >0    will   that      the   locus   of    the    det  (I+G) 
does    not    pass    through  the   -1      point.      If    the   1      or    Euclidean 
norm    is   assumed      for    this    condition    the   equation       3.2    may    be 
expressed   in   terms    of  singular    values  as   equation    3.3. 
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Ov  (AG)*  0\  (I    ♦  G)  (3.3) 

This  result  states  that  as  long  as  the  maximum  singular 
value  of  the  perturbation  matrix  A  G  is  below  the  minimum 
norm  value  of  the  return  difference  matrix  the  system  will 
remain  stable.  The  problem  of  guaranteeing  robustness 
becomes  that  of  finding  the  largest  norm  of  the  perturbation 
guantity,  the  largest  singular  value,  for  which  the  smallest 
norm  or  singular  value  of  the  return  difference  matrix  will 
remain  non-singular. 

The  multiplicative  form  for  a   system  such  as  Figure  3.4 
gives  the  similar  norm  eguation  in  eguation  3.4 

I  |AG(jaJ)i  |<V|  I  (I  +  (GJ-i)-MI  (3.4) 

(JO    >0.  which  may  be  expressed  as 

^  (AG)  <  Qv  (  I  +  G-i  )  (3.5) 

Singular  value  decomposition  software  is  readily  available 
to  determine  how  near  the  matrix  I  +  G_  or  1+  (G)_1  is  to 
singularity. 
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Figure  3.4   Multiplicative  Perturbation. 


25 


IV.  POLE  PLACEMENT  ONLY  DESIGN  PBOCEDOEE 

The  designer  may  use  either  objective,  constraint  or  a 
combination  of  functions  to  secure  the  desired  pole  loca- 
tions. As  currently  implemented  in  the  program  the  cost  or 
objective  portion  of  the  pole  placement  procedure  is 
constructed  as  equation  4.1 

I 

OBJ  =  Z   <  X^-  AR,  )  ^  ♦  (  X  -  \1a  )   2  (4.1) 

A  =1 

where      /\R  =  real  eigenvalue 

Ai  =  imaginary  eigenvalue 

7\rp  =  desired  eigenvalue  location 

7\jo  =  desired  eigenvalue  location 
The  constraint   formulation  is  a   function  that   must  be 
kept  negative  or  the  constraint   is  violated.   It  is  written 
as  equation  4.2 

g(j)  =  ^(Awr^V2  +  i>w-Xi,-)2  -  r  <4-2) 

where  r  is  a  tolerance  circle  established  as  a  function  of 
pole  placement  position.  Since  the  aim  of  the  optimizer  is 
to  keep  g  negative  any  time  the  A  function  of  the 
constraint  is  greater  than  r  the  constraint  will  become 
active,  i.e.  violated.  The  optimizer  will  than  attempt  to 
move  the  constraint  to  the  inactive  status  by  adjusting  the 
design  parameters  of  the  system.  For  obtaining  the  feedback 
gains  and  desired  poles  use  the  'CONXSV  program  [Eef.  9], 
weight  function  1  equal  1  the  other  weight  functions  are  set 
to  0.  After  obtaining  the  feedback  gains  the  NPS  'OPTSYS1 
program  was  used  to  plot  the  Bode,  Nyguist  diagram  and  pole 
zeros  plot. 
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V.  POLE  PLACEMEHT  AHD  ROBUSTNESS  DESIGN  PROCEDURE 

Consideration  of  implementation  of  the  frequency  domain 
or  robustness  portion  of  the  design  procedure  begins  with 
the  concept  of  MIMO  phase  and  gain  margins.  Several  useful 
theorems  on  singular  value  analysis  of  multiloop  systems  are 
presented  in  [Ref.  3]-  One  of  these  theorems  relates  the 
matrix  singular  value  of  the  return  difference  function  to  a 
parameter, $.  ,  and  further  shows  that  as  long  as  the  maximum 
singular  value  of  the  perturbation  function  (L_1  -  I) 
remains  less  than  this  oi.r  the  system  remains  stable.  The 
value  of  o(.  is  then  related  to  gain  and  phase  margins  of  the 
MIMO  system.  The  relationship  developed  is  given  in  equa- 
tions 5. 1  and  5. 2: 

gain  margin  =  GM  =1/(1+o(p)  (5.1) 

phase  margin  =  PM  =  +cos  -»  ( 1-  c(,2/2)  (5.2) 

provided  that  equation  5.3  holds. 

OMI  +  G)>  oi„  (5.3) 

for  some  ck0  <  1 
These  phase  and  gain  margins   are  guaranteed  in   every  loop 
simultaneously. 

Universal  gain  and  phase  margin  curves,  [Ref.  5],  based 
on  the  minimum  singular  values  of  the  return  difference 
matrix  are  developed  from  equation  5.4. 

(L-i-I)  =max  ^  (1-  i/k7)  2  +  2/k„  (1-:os/J  (5.4) 
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foe  all  n   with  k^  >  0. 


These  curves  shown  in  Figure  5.1 
allow  the  designer  to  pick  a  singular  value  that  corresponds 
to  a  specific  gain  and  phase  margin  for  a  given  system. 
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Figure  5.1    Universal  Gain. and  Phase  Singular  Value  Plot. 

Since  the  universal  curve  in  Figure  5.1  provides  a 
convenient  method  of  specifying  gain  and  phase  margins  in 
terms  of  singular  values  the  robustness  portion  of  the  pole 
placement  and  robustness  design  procedure  uses  the  minimum 
singular  value  level  of  the  return  difference  matrix  to 
determine  the  robustness.  The  minimum  singular  value  level 
is  formulated  as  an  objective  or  constraint  function  in 
equation  5.5 


J=  L(max(0,  ((£-  o-(jw,p)))  2 


(5.5) 
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The  optimization  procedure  may  be  used  to  change  feedback 
gains  until  the  minimum  singular  value  is  raised  above  this 
desired  design  level.  Although  the  same  formulation  can  be 
used  as  a  negative  constraint  function  it  has  not  been 
implemented  as  such  within  this  program.  There  are  numerous 
ways  the  singular  value  formulation  could  be  implemented 
vithin  the  program  by  changes  of  the  code  if  design  reguire- 
ments  forced  such  changes. 

The  pole  placement  and  robustness  design  program  is 
based  on  the  ADS  code  to  implement  the  design  variable 
selection  procedures.  The  pole  placement  and  robustness 
program  is  used  to  provide  designs  for  state  or  output  feed- 
back problem. 

For  the  state  or  output  feedback  design  program  the  user 
must  input  the  plant  matrices  A_,B,  C_  and  initial  starting 
values  for  the  feedback  matrix  F.  The  matrices  correspond  to 
the  following  linear  differential  system: 

x  =  Ax  +  Bu  (5.6) 

y  =  Cx  (5.7) 

u  =  -Fx  (5.8) 

As  the  design  program  is  currently  coded  the  user  may 
run  output  feedback  or  state  feedback  design  by  specifying 
the  C  matrix  as  the  diagonal(I)  matrix  for  state  feedback. 
The  program  relies  on  initial  starting  values  of  the  feed- 
back   gains,    F. 

The  ability  to  select  acceptable  starting  values  for  the 
feedback  gains  will  make  the  procedure  more  efficient  in 
operation.  As  currently  employed,  the  program  is  used  to 
obtain  pole  placement  and  robustness  for  a  given  set  of 
starting    gains    and   a    selected      optimization    routine   from    the 
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ADS  program.  If  the  optimizer  is  not  able  to  meet  the 
desired  design  goals  on  this  program  run  two  options  are 
available.  First,  change  to  a  different  optimization  routine 
from  the  list  of  available  ADS  routines  and  rerun  the 
problem.  This  was  usually  successful  in  improving  the 
design.  Second,  the  designer  uses  a  new  set  of  starting 
values  for  the  feedback  gains  and  repeats  the  design  proce- 
dure. Both  options  might  be  used  on  particularly  difficult 
cases. 

The  pole  placement  and  robustness  design  algorithm 
computes  input  additive,  output  adiitive,  input  multiplica- 
tive, and  output  multiplicative  singular  values.  The  versa- 
tility of  the  pole  placement  and  robustness  design  is 
obtained  by  incorporating  a  state  of  the  art  optimizer 
routine  ADS,  with  currently  available  singular  value  compu- 
tation routines. 
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VI.    APPLIC AT  IONS (X22 A    LONGITUDINAL    PROBLEM) 

This  chapter  will  deal  with  a  more  practical  application 
of  the  numerical  optimization  program.  In  this  problem  the 
combined  pole  placement  only,  robustness  design  procedure 
will  be  applied  to  the  linear  longitudinal  dynamic  channels 
of  an  X22A  V/STOL  a/c  [Eef.  8].  The  model  is  that  of  the 
longitudinal  dynamics  of  an  X22A  V/STOL  a/c  at  low  altitude 
and    airspeed   =    65   knots.    The   dynamic    model    of    ttie    system    is 


x    =   A   x    +    B    u 


(6.1) 


X    ={u/wrg,9) 


(6.2) 


£     =    (cC*cf) 


(6.3) 


where 

'  -0.18 
-0.2 
-0.01 
0 


and   where 


B    = 


^ 


-0.36 
0. 

0.33 
0. 


-0.03 
-0.55 
-0.02 
0 

0.52 
-1. 
0.02 
0. 


9.57 
109.43 
-0.  1 
1. 


-31.87 
2.78 
0 
0 


with    full  state    available   for    feedback.    Table    1    is    a    summary 
of    parameters.       The    system    is   not    open-loop   stable. 

The  control  law  is  formulated  to  satisfy  the  desired 
performance  specification.  Equation  6.4  is  the  basic  control 
law. 
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TABLE    1 
X22A    V/STOL  a/c      Parameter   Definitions 


Variable 

Units 

Description 

0 

W 

Q 

e 

ft/sec 
ft/sec 
rad/sec 
rad 
inch 
inch 

Body  axis  forward   velocity  pertui bati> n 
Vertical  velocity   perturbation 
Pitch  rate 
Pitch  attitude 
Elevator  position 
Throttle  position 

u    =   -Pj 


(6.4) 


The      design   studies      presented    up      to  this      point    have      been 
based   on   breaking   the  system   locp  at   the   input   as    shown   in 


Figure   6.1        Systea  Block   Diagram. 
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Figure  6.1.  In  multivariable  tneory  the  location  of  the 
break  in  the  loop  changes  the  return  difference  for  the 
system  and  the  transfer  function  formulation.  In  Figure  6.1, 
number  1  depicts  a  system  with  an  input  loop  break  point 
while  number  2  depicts  an  output  loop  break  point  for  output 
return  difference  determination.  The  return  difference  func- 
tion for  point  1  is  written  as  I+F3  while  the  return  differ- 
ence for  point  2  is  I+GF.  In  figure  number  1  the  transfer 
function  FG  is 

FG  =  F'*  (SI  -  A)  »*E 

F*  =  F*C 

and  GF  is    GF  =  C*(SI  -  A)   **B» 

Bf  =  B*F 


A.       INPUT   TEANSFEB    FUNCTION    F*G   ANALYSIS 

Using  the  NPGS  'CONXSV*  program  [Eef.  9],  to  get  the 
desired  feedback  constants,  the  pole  placement  only  program 
was      first   run      to   get     the      desired   pole      and   the      computed 


TABLE   2 
Pole   Only    Desired    Pole  and   Computed  Pole 


Desired   pole 

Computed   pole 

-1.2               1.61 
-1.2            -1.63 

-1.0                      J 
-0.5 

-1.19869                   1.594151 
-1.19889                -1. 59415 j 
-0.98902 
-0.49770 

pole  (see  Table    2)  .       The   feedback   matrix    F    is 

-0.10853       0.50675       13.97235         7.73604 
0.19809    -0.10997   -90.41902      12.42874 
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Figure  6.2    Pole  Placement  Only  Singular  Value  Plots. 


the  objective  function  -  is  0.0002  and  the  additive  minimum 
singular  value  vs  frequency  is  seen  in  Figure  6.2.  This 
plot  shows  that  pole  placement  only  design  has  very  low 
minimum  singular  values  for  the  return  difference  matrix. 
Pole  placement  only  goes  as  low  as  -22  db  or  0.0827  rad  near 
2  rad/sec  in  frequency.  Using  the  universal  gain  and  phase 
diagram  as  discussed  in  Chapter  5  this  equates  to  a  gain 
margin  of  about  0.95  db  to  1.15  db  and  a  phase  margin  of  4.7 
degrees.  These  phase  and  gain  margins  are  quite  small, 
showing  the  need  to  run  a  pole  placement  and  robustness 
design  program.  It  is  assumed  that  the  eigenvalues  (pole 
locations)  of  the  system,  as  developed  in  [Bef.  7],   are  the 
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required  poles  for  the  performance  criteria.  Once  the  pole 
locations  are  set,  robustness  criteria  must  be  selected. 
From  the  universal  gain  and  phase  margin  curve  discussed 
earlier,  singular  value  levels  were  selected  for  this 
problem.  The  singular  value  level  chosen  was  0.6  rad.  This 
corresponds  to  a  gain  margin  of  -4.0  db  to  8  db  and  a  phase 
margin  of  about  35  degrees. 

For  a  singular  value  level  cf  0.6  rad  the  pole  placement 
and  robustness  design  routine  places  the  poles  as  shown  in 
Table  3.     The  slight  differences   in  these   pole  locations 

TABLE  3 
Pole  Placement  and  Bobustness  Design  Pole  Locations 


Desired   pole 

Computed    pole 

-1.2               1.61 
-1.2             -1.6J 
-1.0 
-0.5 

-1.19297                   1.595331 
-1.  19297               -1.595331 
-0.99886 
-0.45480 

have  an  insignificant  effect  on  the  performance. 
The  feedback  constant  F  is 


F   = 


-0.07493   -0.08746    8.0529 
0.27085   -0.23019   -2.2919 


16.30763 
37.66504 


/ 


The  objective  function  is  0.0000  and  the  additive 
minimum  singular  value  vs  frequency  is  seen  in  Figure  6.3. 
This  plot  shows  that  in  the  robustness  design  the  gain 
adjustment  moves  the  minimum  singular  value  from  about  0.08 
rad,  with  very  poor  phase  and  gain  margins,  to  a  level  of 
0.645   rad,  above    the      desired   values      of   gain      and    phase. 

After  obtaining  the  feedback  gains  the  NPGS  'OPTSYS'  program 
was  used  to  obtain  the  necessary  data,  time  response  plot, 
Eode  plot,  and  pole-zero  map.  Ths  pole  placement  only 
design   closed      loop    time      response    plot      is    shown      in   Figure 
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Figure  6.3    Robustness  Design  Singular  Value  Plots. 

6.4.  It  shows  good  pitch  rate  and  pitch  attitude  responses 
to  (£  step  input,  but  slow  responses  to  the  forward 
velocity  perturbation  and  vertical  velocity  perturbation  to 
Ot  step  input.  The  input  transfer  function  F*G  robustness 
design  time  response  plot  is  shown  in  Figure  6.5.  This 
indicates  a  better  time  response  for  pitch  rate  and  pitch 
attitude  inputs  than  pole  placement  only  design,  also  good 
vertical  velocity  perturbation  response  to  oQ  step  input, 
but  sluggish  response  for  the  forward  velocity  perturbation 
[Bef.  8].  Comparing  the  loop  Eode  plots  of  the  pole  place- 
ment and  robustness  design  transfer  function,  F*G,  shows  an 
increase  in  robustness. 
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Figure   6.4        Pole   Placeaent    Only   Design    Time    Response. 
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Figure  6.5    F*G  Robustness  Design  Time  Response. 
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Considering  a  SISC  Bode  analysis,  the  channel  input  1 
output  1  bode  plot  are  shown  in  Figure  6.6  and  6.7. 
The  bandwidth (BR)  shows  a  small  decrease,  3  rad/sec,  and  an 
increase,  in  the  phase  margin (EM),  of  22  degrees.  The  DC 
gain  drops  from  34.3  db  to  25.7  db,  a  9  db  decrease  and  the 
slope  changes  very  abruptly  in  the  pole  only  design  (see 
Figure  6.4).  The  robustness  design  seen  in  figure  6.5  drops 
the  curve  down  20  db  due  to  the  two  zeros  at  0.3  rad/sec, 
which  move  toward  the  two  poles  at  1.27  rad/sec,  and  the 
other  zero  moves  closer  to  another  pole  0.16  rad/sec  during 
the  optimization.  This  channel  only  shows  a  smoothing  of  the 
curve.  Consider  the  channel  input  1  output  2  bode  plots 
shownin  Figures  6.8  and  6.9.  This  channel  shows  a  large 
shift  in  the  bandwidth (EU)  from  30  rad/sec  to  2.83  rad/sec. 
The  DC  gains  are  similar  and  increase  the  gain  margin  24  db, 
but  decrease  phase  margin  90  degrees.  Figure  6.3  shows  that 
the  frequency  response  curve  slope  changes  detrimentally, 
but,  the  robustness  design,  as  shown  in  figure  6.9  shows  a 
decrease  in  the  slope  down  20  db  curve  by  zeros  move  the  one 
pole  0.17  rad/sec  for  pole-zero  cancel  and  the  other  one 
zero  moves  on  the  minimum  singular  value  frequency  1 
rad/sec.  The  reduction  of  bandwidth  and  increase  in  gain 
margin  yields  an  increased  tolerance  to  perturbation.  The 
channel  input  2  output  1  bode  plots  are  shown  in  Figures 
6.10  and  6.11.  The  DC  gain  decreases  10  db  and  bandwidth 
shifts  left  0.1  rad/sec,  increasing  the  phase  margin  43.7 
degrees.  The  zero  shift  has  the  effect  of  smoothing  the 
frequency  response  curve  in  the  vicinity  of  the  frequency  of 
the  minimum  singular  value.  The  channel  input  2  output  2 
bode  plots  are  shown  in  Figures  6.12  and  6.13.  These 
figures  show  the  same  response  as  the  channel  1;2  but  a 
bandwidth  shift  right  of  0.47  rad/sec  and  the  same  DC  gainas 
for  a  smooth  curve,  caused  by  the  zero  location  shift  in  the 
same   manner  as   the   other  channels (channels   1 ; 1 , 2; 1 , 2; 2) . 
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Figure   6.6        Pole   Placement   Only    Design      F*G    1:1, 
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Figure    6-7        Robustness   Design      F*G    1:1 
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Figure  6.8        Pole  Placement  Only   Design      F*G    1:2. 
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Figure   6.9        Robustness   Design      F*G    1:2. 
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The  r*G  transfer   functions  bode  plot  numerical   results  are 
shown  in  Table  4. 

TABLE  4 
F*G  Transfer  Functions  Bode  Plot  Results 


Channel 

Bandvidth 
rad/sec 
pole     SV 

PM 
degree 
pole   SV 

GM 

Db 

pole 

SV 

Gain 
Db 
pole 

SV 

1  : 

:  1 

6.1     3.2 

58    80 

• 

• 

34 

25.7 

1  : 

:  2 

30      2.83 

89.7  -2  0 

-32 

-7.2 

28 

27.6 

2  : 

:  1 

0.3     0.19 

9.5   53 

8.8 

• 

19.4 

8.6 

2 

:  2 

0.  17    0.64 

170   106 

-7.3 

• 

13.5 

14.1 

The  most  noticed  change  in  the  overall  system,  however, 
is  in  the  transfer  function  input  1  output  2  .  This  is  the 
change  at  the  minimum  singular  value  position,  which  greatly 
reduces  the  gain  and  bandwidth  in  this  channel  through  a 
change  in  feedback  gain.  The  optimizer  routine  Drings  the 
entire  system  gains  to  more  balanced  conditions  and  recovers 
a  highly  robust  design. 

The  gain  changes"  associated  with  the  robustness  improve- 
ment cause  the  zeros  of  the  various  closed-loop  pole-zero 
diagram  of  the  closed-loop  transfer  matrices  to  move.  A 
comparison  of  the  eight  pole-zero  diagrams  is  shown  in 
Figures  6.14  to  6.17. 

The  significant  feature  of  these  pole-zero  diagrams  is 
the  shift  in  the  zeros  of  the  optimized  design  in  a  direc- 
tion that  attempts  to  equalize  or  balance  the  frequency 
response  for  frequencies  in  the  vicinity  of  the  minimum 
singular  values.  The  pole-zero  diagram  of  channel  input  1 
output  2  will  be  discussed  as  an  example  of  this  effect.  In 
figure  6.12  the  pole  only  design  zeros  are  located  at  about 
-0.31  and  -0.78.  When  the  pole  placement  and  robustness 
routine  has   completed  the  feedback  gain   modification  these 
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zeros  have  shifted  to  -0.17  and  -0.56.  The  effect  of  these 
zero  shifts  is  to  combine  with  the  pole  locations  to 
equalize  the  frequency  response  and  increase  DC  gain  of  14 
db  for  the  same  bandwidth,  at  the  minimum  singular  value 
freguency  1  rad/sec  as  depicted  in  Figures  6.18  and  6.19. 
Zero  shifts  for  the  remainder  of  the  transfer  functions 
provide  similar  results  in  the  other  channels.  By  moving 
toward  the  frequencies  associated  with  the  minimum  singular 
values  the  zeros  have  balanced  the  overall  freguency 
response  of  the  system  in  each  channel.  While  the  channel 
gain  modification  is  the  primary  mechanism  for  robustness 
recovery,  the  zero  shift  associated  with  the  feedback  gain 
changes  is  directly  related  to  the  overall  frequency 
response  of  the  system.  A  robustness  design  meeting  the 
reguired  pole  locations  and  a  robustness  singular  value 
level  of  0.6  rad  provides  adeguate  gain  and  phase  margin  for 
the  design. 


B.   ODTPDT  TRANSFER  FUNCTION  G*F  ANALYSIS 

Using  the  NPGS  'CONXSV1  program   [Eef.  9],   to   get  the 

desired  output  transfer  function   feedback  constants,    the 

pole  placement  only  program  was  first  run  to  get  the  desired 

pole   and  the  computed  pole (see  Table   2).    the   feedback 

matrix  F  is 

s 

-0.10853       0.50675       13.97235         7.73604 

0.19809    -0.10997    -90.41902       12.42874 

the    objective      function  is   0.0002  and    the      additive   output 

minimum    singular    value  vs    frequency  is   seen    in   Figure   6.20. 

This    plot   shows    that    pole   placement  only   G*F   dssign    has   very 

low      minimum      singular      values      for  the      return      difference 
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Figure   6.18        Pole   Placement    Only   Design    F*G   Closedloop    1:2. 
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Figure  6.19    Robustness  Design   F*G  Closedloop  1:2 
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Figure  6,20    Pole  placenent  only  3*F  Singular  Value  Plots. 

matrix.  Pole  placement  only  goes  as  low  as  -48  db  or  0.004 
rad  near  1  rad/sec  in  frequency.  Using  the  universal  gain 
and  phase  diagram  as  discussed  in  Chapter  5  this  equates  to 
a  gain  margin  of  about  -0.003  db  to  0.003  db  and  a  phase 
margin  of  0.5  degrees.  These  phase  and  gain  margins  are 
quite  small,  showing  the  need  to  run  a  pole  placement  and 
robustness  design  program.  It  is  assumed  that  the  eigenva- 
lues (pole  locations)  of  the  system,  as  developed  in 
[Ref.  7],  are  the  required  poles  for  the  performance 
criteria.  Once  the  pole  locations  are  set (see  Table  5), 
robustness  criteria  curve  discussed  earlier,  singular  value 
levels  are  made  for  this  problem.  The  singular  value  level 
chosen  is  0.6  rad.  This  corresponds  to  a  gain  margin  of  -4.0 
db  to  8  db  and  a  phase  margin  of  about  35  degrees. 
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For  a  singular  value  level  of  0.6  rad  the  pole  placement 
and  robustness  design  routine  places  the  poles  as  shown  in 
Table  5.  There  is  a  trade  off  between  performance  and 
robustness  in  many  control  problems.  IN  output  transfer 
function  G*F  singular  value  calculation,  the  shift  in  this 
pole  can  have  a  significant  effect  on  the  performance,  it 
was  found  necessary  to  move  the  pole  at  -0.5  to  -0.05  in 
order  to  achieve  the  desired  rotustness.  The  slight  differ- 
ences in  these  pole  locations  an  insignificant  effect  on  the 

TABLE  5 
G*F  Pole  Placeaent  and  Robustness  Design  Pole  Locations 


Desired   pole 

Computed    pole 

-1.2               1.61 
-1.2             -1.6 i 
-1.0 
-0.05 

-1.20944                    1.58904J 
-1.20944                -1.589041 
-1.02955 
-0.05836 

performance.   The  feedback  matrix  F  is 


/ 


-0.27091         0.18869         0.0 
1.88298      -1.60017        0.0 
The     objective      function      is  0.0002 


0.0 

0.00914 

and   the   additive 


output  singular  value  vs  frequency  is  seen  in  Figure  6.21. 
This  plot  shows  that  in  the  robustness  design  the  gain 
adjustment  moves  the  minimum  singular  value  from  about  0.004 
rad,  with  very  poor  phase  and  gain  margins,  to  a  level  of 
0.612  rac,  close  enough  to  the  desired  values  of  gain  and 
phase.  After  obtaining  the  feedback  gains  the  NPGS  'OPTSYS' 
program  was  used  to  obtain  the  necessary  data,  Bode  plot, 
pole-zero  map,  and  time  response  plot.  The  G*F  robustness 
design  time  response  plot  is  shown  in  Figure  6.22.  The 
response  of  pitch  rate  is  relatively  good  but  sluggish 
response  of  pitch  attitude  to  c£  step  input.  The  performance 
of   the    pitch      attitude   was    degraded.      It    is      possible    that   a 
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Figure  6.21    Robustness  Design  G*F  Singular  Value  Plots. 

further  change  in  the  pole  location  could  achieve  better 
performance  and  still  meet  the  robustness  criteria.  The 
body  axis  forward  velocity  and  vertical  velocity  perturba- 
tion responses  are  the  same  as  in  the  pole  placement  only 
design(see  Figure  6.4).  Comparing  the  loop  Bode  plots  of 
the  pcle  placement  and  robustness  design  transfer  function, 
G*F,  shows  an  increase  in  robustness. 

Considering  a  SISO  Bode  analysis,  the  channel  input  1 
output  1  bode  plot  are  shown  in  Figure  6.23  and  6.24.  The 
bandwidth (BW)  shows  a  small  increase,  0.6  rad/sec,  and  the 
phase  margin  decrease  of  26  degrees.  The  DC  gain  drops  from 
21.4  db  to  17.3  db,  a  4  db  decrease  and  the  slope  change  is 
similar  to  the  pole  only  design  as  seen  in  Figures  6.23  and 
6.24.   This  channel  only  shows  a  smooth  curve. 
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Figure  6.22    G*F  Robustness  Design   Time  Besponse. 
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Figure   6.23         Pole   Placement    Only   Design      G*F    1:1 
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Figure  6.24    Robustness  Design   G*F  1:1 
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Considering  the  channel  input  1  output  2  with  bode  plots 
shown  in  Figures  6.25  and  6.26,  a  shift  in  the  bandwidth  (3W) 
from  2.4  rad/sec  to  3.2  rad/sec  is  shown.  The  DC  gain 
increases  from  12.7  db  to  15.8  db  at  the  minimum  singular 
frequency  of  1  rad/sec  and  an  increase  in  the  phase  margin 
of  28  degrees.  Figures  6.25  and  6.26  show  the  similar 
frequency  response  curve,  slope  changes  by  zeros  move  the 
one  pole  position  0.8  rad/sec  for  pole-zero  cancel  and  the 
other  one  zero  move  to  close  the  minimum  singular  value 
frequency  1  rad/sec  and  two  pole  location.  The  increase  in 
the  DC  gain  and  phase  margin  yields  an  increased  tolerance 
to  perturbation.  The  channel  input  2  output  1  bode  plots 
are  shown  in  Figures  6.27  and  6.28.  The  DC  gain  decreases 
16  db  and  the  bandwidth  shifts  left  0.8  rad/sec,  increasing 
the  phase  margin  151  degrees.  The  zero  shift  has  the  effect 
of  smoothing  the  frequency  response  curve  in  the  vicinity  of 
the  frequency  of  the  minimum  singular  value.  The  channel 
input  2  output  2  bode  plots  are  shown  in  Figures  6.29  and 
6.30.  These  figures  show  a  small  shift  left  in  the  band- 
width from  4.5  rad/sec  to  2.8  rad/sec.  The  DC  gain  decreases 
12  db  with  an  increased  phase  margin  of  56  degrees.  This 
results  in  a  smooth  curve,  caused  by  the  zero  location  shift 
in  the  same  manner  as  previous  channels (channels 
1; 1# 2; 1, 2;2) .  The  G*F  transfer  functions  bode  plot  numer- 
ical results  are  provided  in  Tatle  6. 

The  noticeable  change  in  the  overall  system,  however,  is 
in  the  transfer  function  input  1  output  2  .  This  is  the 
change  at  the  minimum  singular  value  position,  increasing 
the  gain  and  phase  margin  through  a  change  in  feedback  gain. 
The  optimizer  routine  brings  the  system  gains  to  a  more 
balanced  condition  and  recovers  a  highly  robust  design. 

The  gain  changes  associated  with  the  robustness  improve- 
ment cause  the  zeros  of  the  various  closed-loop  pole-zero 
diagram   of  the   closed-loop   transfer   matrix  to   move.    A 
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Figure  6.25    Pole  Placement  Only  Design   G*F  1:2. 
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Figure  6.26   Robustness  Design   G*F  1:2, 
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Figure   6.27        Pole   Placement    Only    Design      G*F    2:1 
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Figure  6.30    Robustness  Design   G*F  2:2, 
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TABLE  6 
G*F  Transfer  Functions  Bode  Plot  Results 


channel 

Bandwidth 
rad/sec 
pole     SV 

PM 
degree 
pole   SV 

GH 
Db 
pole    SV 

Gain 
Db 
pole    SV 

1  :  1 

1  :  2 

2  :  1 
2  :  2 

0.95    1.56 
2.36    3.2 
2.2     1.43 
4.5     2.8 

49    23 
211.3  239 
80   231 
10.7  17.5 

3.9 

•  • 

•  • 

•  • 

21.4   17.3 
12.6   15.8 
37.3   21.5 
29.1   17.5 

comparison   of  the   eight   pole-zero   diagrams  is   shown  in 
Figures  6.31  to  6.34. 

The  significant  feature  of  these  pole-zero  diagrams  is 
the  shift  of  the  zeros  of  the  optimized  design  in  a  direc- 
tion that  attempts  to. equalize  or  balance  the  frequency 
response  for  frequencies  in  the  vicinity  of  the  minimum 
singular  values.  The  pole-zero  diagram  of  channel  input  1 
output  2  will  be  discussed  as  an  example  of  this  effect.  In 
figure  6.30,  the  pcle  only  design  zeros  are  located  about 
-0.31  and  -0.78.  When  the  pole  placement  and  robustness 
routine  completes  the  feedback  c,aiu  modification  these  zeros 
shift  close  to  -1.7  and  0.28.  The  effect  of  these  zeros 
shifting  is  to  combine  with  the  pole  locations  to  equalize 
the  frequency  response  and  increase  the  DC  gain  25  db  at  the 
same  bandwidth  and  the  minimum  singular  value  frequency  1 
rad/sec,  as  depicted  in  Figures  6.35  and  6.36.  Zero  shifts 
for  the  remainder  of  the  transfer  functions  provide  similar 
results  in  the  other  channels.  By  moving  toward  the 
frequencies  associated  with  the  minimum  singular  values  the 
zeros  have  balanced  the  overall  frequency  response  of  the 
system  in  each  channel.  While  the  channel  gain  modification 
is  the  primary  mechanism  for  robustness  recovery  the  zero 
shift  associated  with  the  feedback  gain  changes  is  directly 
related  to  the  overall  fre3uency  response  of  the  system. 
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Figure   6.31        Pole-Zero   Map   for    1:1    and    1:2 
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Figure  6.32    Pole-Zero  Map  for  1:3  and  1:4. 
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Figure    6.33        Pole-Zero   Map    for    2:1    and    2:2, 
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Figure  6.34    Pole-Zero  Map  for  2:3  and  2:4, 
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Figure   6.35        Pole    Placement    Only   Design    G*F   Closedloop    1:2. 
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Figure  6.36    Robustness  Design   G*F  Closedloop  1:2 
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To  summarize  this  analysis,  a  given  performance  level 
has  been  chosen  in  terms  of  pole  locations.  The  level  of 
robustness  has  been  set  foe  a  desired  gain  and  phase  margin 
based  on  the  universal  gain  and  phase  margin  curve.  The  pole 
placement  and  robustness  routine  improves  the  robustness 
level  by  changing  the  feedback  gains  that  affect  the  channel 
input  1  output  2  cross-coupling.  This  robustness  recovery 
is  affected  by  modification  of  the  system  feedback  gains  in 
such  a  manner  that  cross  coupling  gains  are  reduced  so  that 
small  cross-coupling  perturbations  do  not  drive  the  system 
into  instability.  The  open-loop  transfer  function  plots  are 
used  to  indicate  how  this  mechanism  operates  and  have  been 
shown  to  be  an  alternative  indicator  of  channels  that  may  be 
affected  by  crossfeed  perturbations.  The  pole-zero  diagrams 
of  the  closed-loop  transfer  functions  of  the  transfer  matrix 
further  indicate  that  zero  movement  is  in  a  direction  that 
equalizes  the  gain  level  of  the  frequency  response  curves  in 
the  vicinity  of  the  lowest  singular  values  providing  a  more 
balanced  system  response.  Finally,  a  comparison  of  the 
input  transfer  function  F*3  output  singular  value  and  output 
transfer  function  G*F  input  singular  value  plots  are  shown 
in  Figures  6.37  and  6.38.  The  F*G  output  minimum  singular 
value  is  too  low(  0.0076  rad).  The  G*F  input  minimum 
singular  value  is  also  its  own  low  (0.1  rad).  Therefore 
both  transfer  functions  control  only  own  function's  robust- 
ness. Further  analysis  would  be  needed  if  robustness  wtre 
required  in  both  the  input  and  output  case  simultaneously. 
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VII.  CONCIOSIONS 

The  pole  placement  and  robustness  design  routine  coupled 
with  the  Automated  Design  Synthesis  program  provides  the 
designer  an  excellent  tool  with  which  to  attack  the  robust 
design  problem.  The  pole  placement  and  robustness  design 
routine  has  demonstrated  the  capability  of  providing  designs 
that  solve  the  X22A  V/STOL  A/C  longitudinal  dynamic  problem 
caused  by  cross-coupling  perturbations  which  reduce  robust- 
ness in  multivariable  systems.  This  design  improvement  is 
accomplished  by  modifying  the  system  feedback  gains  in  such 
a  manner  that  the  gain  in  channels  that  are  affected  by 
cross-coupling  perturbations  is  equalized  with  other  system 
gains  to  reduce  this  cross-coupling  effect.  The  gain  changes 
are  accompanied  by  zero  shifts  which  also  influence  the  gain 
distribution  and  freguency  response  of  the  system. 

The  lack  of  robustness  can  be  discovered  in  two  ways. 
The  first  method  is  to  plot  the  open-loop  Bode  plots  of  each 
element  of  the  transfer  matrix  and  look  for  extremely  high 
gains  or  bandwidths  relative  to  the  other  transfer  func- 
tions. The  second  method  examines  the  singular  values  of 
the  return  difference  matrix  for  magnitude.  The  joint  anal- 
ysis method  tells  where  the  robustness  problem  occurs  in  the 
design.  Low  singular  values  correspond  to  low  robustness. 
The  pole  placement  and  robustness  design  routine  can 
increase  robustness  by  modifying  feedback  gains  to  reduce 
the  effect  of  cross-coupling  within  the  system.  Observing 
the  gain  modification  made  by  the  pole  placement  and  robust- 
ness routine  the  critical  channel  within  the  system  that 
affects  the  robustness  can  be  determined  from  the  Bode 
plots.  The  pole  placement  and  robustness  routine  feedback 
gain  changes  also  cause  zero   shifts  during   the  robustness 
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recovery.  The  gain  on  the  open  loop  Bode  plot  for  the 
affected  cross-coupling  channel  is  adjusted  and  the  closed- 
loop  zeros,  as  seen  on  the  pole-zero  diagram,  are  shifted. 
This  zero  shift  is  in  a  direction  which  will  combine  with 
system  poles  to  smooth  the  frequency  response  jliwigram  in  the 
vicinity  of  the  minimum  singular  value.  The  performance  of 
the  robustness  design  input  transfer  function  F*G  provides 
the  best  response,  but  the  output  transfer  function  G*F 
indicates  the  performance  of  the  pitch  attitude  is  degraded. 
It  is  possible  that  a  further  change  in  the  pole  location 
would  provide  better  performance  and  still  meet  the  robust- 
ness   criteria. 
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